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Abstract 

The effective thermal conductivity of high-porosity 
open cell nickel foam samples was measured over a 
wide range of temperatur es and pressures using a 
standar d steady-state technique. The samples, 
measuring 23.8 mm, 18.7 mm, and 13.6 mm in 
thickness, were constructed with layers of 1.7 mm thick 
foam with a porosity of 0.968. Tests were conducted 
with the specimens subjected to temperature differences 
of 100 to 1000 K across die dtickness and at 
environmental pressures of 10 4 to 750 mm Hg. All test 
were conducted in a gaseous nitrogen environment. A 
one-dimensional finite volume numerical model was 
developed to model combined radiation/conduction 
heat transfer in the foam. The radiation heat transfer 
was modeled using the two-flux approximation. Solid 
and gas conduction were modeled using standard 
techniques for high porosity media. A parameter 
estimation technique was used in conjunction with the 
measured and predicted thermal conductivities at 
pressures of 10 4 and 750 mm Hg to determine the 
extinction coefficient, albedo of scattering, and 
weighting factors for modeling the conduction thermal 
conductivity. The measured and predicted 
conductivities over the intermediate pressure values 
differed by 13 percent. 


Introduction 

Thermal Protection Systems (TPS ) are currently under 
development for a wide range of aerospace 
applications. Metallic foam has been considered as the 
insulating material of the TPS on reusable launch 
vehicles. Another application is to use the metallic 


foam as part of an integrated structure that serves as 
the launch vehicle's primary structure and thermal 
protection system. However, limited information is 
available on the thermal properties of metallic foams 
under the environmental conditions to which re-entry 
type vehicles are exposed. Earth re-entry typically 
produces aerodynamic heating to a surface to 
temperatures as high as 1000°C in a pressure range 
from 10 2 to 760 mm Hg. 1 

Metallic foams have been recently investigated for 
various applications. Aluminum foams have been 
utilized in the automobile industry for their crash worthy 
properties. 2 Nickel foams have been used to improve 
the performance in high-power batteries. Foam-based 
nickel metal hydride batteries currently compete with 
the more expensive lithium ion batteries for lightweight 
cordless electronics.’ 

There has been extensive work in both experimental 
and analytical modeling of heat transfer in porous 
media. Kaviany has provided a comprehensive review 
of heat transfer in general porous material 4 Lee and 
Cunnington have provided an extensive review of 
conduction and radiation heat transfer in high porosity 
fibrous insulation. 5 Glicksman has reviewed heat 
transfer in polymeric foams. 1 ’ Gibson and Ashby 
discussed thermal properties of foams in their 
comprehensive work on cellular solids/ Ballis, et. al., 
modeled heat transfer in open cell car bon foams and 
determined radiation as the primary source of heat 
transfer for temperatures above 1000 K. 8 They used the 
optically thick approximation for modeling radiation 
with a weighted spectral extinction coefficient to 
account for anisotropic scattering, hi addition, they 
used lineal' superposition of solid conduction, gas 
conduction, and radiation thermal conductivities. Their 
work was limited to atmospheric pressure. There has 
been limited work on metallic foams. Calmidi and 
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Maliajan measured the thermal conductivity of 
aluminum foams at temperatures up to 75° C with a 
15°C temperature difference maintained across the 
specimen. 9 Under these experimental conditions 
nonlinear effects such as natural convection and 
radiation were determined to be negligible. Their work 
was also limited to atmospheric pressure. They modeled 
die heat transfer as combined gas and solid conduction 
and formulated an empirical correlation for die thermal 
conductivity in die aluminum foam. Bhattachrya, et. 
al., performed further refinements on die gas/solid 
conduction modeling of Calmidi and Maliajan, again 
ignoring radiation heat transfer. 10 They showed diat 
die effective thermal conductivity had a strong 
dependence on die porosity and die particular 
geometrical configurations of the intersection of the 
struts hi die foam. 

The objective of the present study was to investigate 
heat transfer hi metallic foams over a wide range of 
pressures and temperatures where solid conduction, gas 
conduction, and radiation were die three modes of heat 
transfer. Therefore, die effective thermal conductivity 
of die nickel foam was measured hi die temperature 
range of 300 to 1300 K and environmental pressures of 
10 4 to 750 mm Hg. The two-flux approximation was 
used to model die radiation heat transfer through die 
nickel foam. A combined conduction formulation 
based on the superposition of solid conduction and gas 
conduction was used to model conduction heat transfer. 
The temperature dependent gas conduction model was 
applicable over the rarefied, transition, and continuum 
gas transport regimes. Inverse heat transfer methods 
were used to determine the parameters needed hi the 
heat transfer model. 


Metallic Nickel Foam 

The metallic nickel foam used hi this study is 
commercially available. It is manufactured hi bulk for 
use hi a variety of applications, hi the manufacturing 
process, polyurethane foam is used as a template. A 
proprietary chemical vapor decomposition process 
coats the surface of the template with nickel. The 
material is annealed at around 1800°C, caushig 
evaporation of the polyurethane core. Photomicrograph 
images of the foam produced by an electron scanning 
microscope with magnification factors of 75 and 750 
are shown in Figures 1 and 2, respectively. As seen 
from the images, the foam has an open cell structure 
with randomly oriented struts, hi addition, the struts 
are themselves hollow, left void due to the removal of 
the polyurethane template. The metallic foam used in 
this study is 99.98% nickel by 



Figure 1. Photomicrograph of nickel foam 
(magnification factor = 75). 


composition with porosity of 0.968. The foam is 1.7 
mm thick and has a density of 290 kg/m 1 at room 
temperature and atmospheric pressure. Although it 
would have been ideal to have samples of different 
thicknesses, the 1.7 mm thick foam was the only 
material available that provided the desired porosity 
and composition. As a result, three samples were 
constructed utilizing 14, 1 1, and 8 layers of the metallic 
foam. The corresponding thicknesses of the samples 
were 23.8 mm, 18.7 mm, and 13.6 mm, respectively. 



Figure 2. Photomicrograph of nickel foam 
(magnification factor = 750). 


Experimental Procedure 

The apparatus used to determine the effective thermal 
conductivity of the foam has been used in previous 
work to investigate heat transfer through fibrous 
insulation. 11 The sample was placed between a 
radiantly heated septum plate that can reach 
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temperatures as high as 1300 K and a water-cooled 
plate maintained at room temperature. The water- 
cooled plate was placed on the bottom so the direction 
of the applied heat flux was opposite to the local gravity 
vector thus eliminating natural convection. 
Thermocouples were placed within die septum plate 
and widiin die water-cooled plate to measure die plate 
temperatures. Heat flux gauges, located on die water- 
cooled plate, measured die flux of heat energy flowing 
dirough die sample. The entire apparatus was housed in 
a 1.5x1. 5 meter vacuum chamber. Upon insertion of 
die foam sample into die test apparatus, shown in 
Figure 3, die gas inside die vacuum chamber was 
removed. Once high vacuum was achieved, nitrogen 
gas was used to regulate die gas pressure widiin die 
chamber. The primary reason nitrogen gas was selected 
for use in diis experiment was to eliminate die 
oxidation of die foam at higher exposure temperatures. 



die 304.8 mm square sample. Widiin die inner region, 
measurements were taken at four different locations. 
Values were recorded at a 5 -second interval for two 
minutes. The average of the measurements of the four 
spatial locations over die two-minute interval was taken 
to be the steady-state measurement. Steady-state 
measurements were taken at 7 different nominal septum 
plate temperatures (from 100°C to 1000°C) and 10 
different gas pressures (from 10 3 4 to 750 mm Hg). 


Analytical Model 

Existing techniques were utilized, modified, and 
combined to model die effective thermal conductivity 
of the metallic foam in die temperature and pressure 
range covered by the experimental data. 1,4 " 11,13 

The conservation of energy principle for a control 
volume, states diat die time rate of increase in die 
energy stored in die volume plus die net rate at which 
energy is conducted out through die surface of the 
volume is equal to zero. The statement can be written 
in integral form as 

Slip- ^■‘ IR+ § 1 " ds = 0 ra 

R Ul S 

where p is die density, c is the specific heat of the 
material, T is the absolute temperature, q is the heat 
flux vector, and n is the normal to the control volume 
surface. 17 


Figure 3. Schematic of experimental apparatus. 

Each measurement consisted of setting the septum plate 
to die desired temperature, varying the nitrogen gas to 
the desired pressure, and allowing both to reach a 
steady state condition. Additionally, the time rate of 
change of the heat flux through the sample was 
monitored until it approached zero. Once steady-state 
environmental conditions were achieved, the 
temperatures of the septum plate and water-cooled plate 
and the heat fluxes were recorded. The effective 
thermal conductivity of the metallic foam, k e jf, was 
determined using Fourier’s law 

k if = — — — (D 

e JJ 'T< rp 

1 Hot 1 Cold 

where L is the thickness of the sample, T Hol is the 
temperature on the hot side of the sample, T Co id is the 
temperature on the cold side of the sample, and q is the 
steady-state heat flux through the sample. Lateral heat 
losses due to edge effects were neglected by restricting 
measurements to the inner 101.6 nun square region of 


When the heat flux vector is unidirectional and parallel 
to the normal of the surface of the control volume (heat 
flow in one direction only), Equation 2 reduces to 

pc ^~fr +q °“' =q '" (3) 

The heat flux into the volume, q, n , and the heat flux out 
of the volume, q out , take the general form 

C l= ( l C oM+ C lra d ( 4 ) 

where q con j is the amount of heat transferred by thermal 
conduction and q rac i is the amount of heat transferred by 
radiation. 

Analytical models for the heat transferred by 
conduction and radiation are discussed and quantified. 
Tire one-dimensional finite volume numerical method is 
presented. Tire nonlinear parameter estimation method 
used to determine intrinsic material properties needed 
in tire analytical formulation is also briefly described. 
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Conduction 


The conduction of heat energy through the foam is a 
combination of conduction through the solid shuts of 
die foam and conduction through die gas widiin die 
void of the foam. The heat flow due to conduction is 
given by 

( 5 ) 

where k is die combined conduction thermal 
conductivity and x is die hi dependent spatial variable in 
die direction of die sample diickness. Several dieories 
have been developed to describe die combined 
conduction diermal conductivity in terms of die 
conductivity of the solid material diat makes up die 
foam and die conductivity of the gas diat fills die voids. 
The limiting case, a parallel arrangement, is based on 
die fraction of cross sectional area responsible for die 
conduction of heat energy. 1 ’ The porosity, e, defined as 
die total void volume divided by die total volume of the 
foam, is used as die fraction of die area diat will 
conduct heat due to gas conduction. Then, (1-f) is die 
fraction of die area that will conduct heat due to solid 
conduction. The parallel arrangement for die combined 
conduction diermal conductivity is given by 

k P ara, M = £ - k g a S +^- £ )- k soM ( 6 ) 


where k gas is the gas diermal conductivity and k so/ld is 
the solid diermal conductivity. Gibson and Ashby 
assumed that the solid thermal conductivity is related to 
the thermal conductivity of the strut’s parent material 
through 


■ F-k, 


(7) 


where F is the solid conduction efficiency factor which 
allows for the tortuous path for conduction through the 
cell walls. 7 The temperature dependent thermal 
conductivity of bulk nickel, kb„ik, is used in Equation 
7. 14 To correct for overestimation of Equations 6 and 
7, Calmidi and Mahajan used the following relationship 
for modeling heat transfer through aluminum foam 

*- = *•*.. + (1 -eY-D-k* ( 8 ) 


and they found that D=0. 18 1 and m= 0.763 produced the 
best match for their experimental results. 9 When D= 1 
and rn= 3, Equation 8 takes the same functional form 
used by Daryabeigi for the combined conduction 
thermal conductivity in fibrous insulation. 11 


Other models investigated assume the conduction takes 
place in a combined parallel-series arrangement of the 
solid and gaseous constituents. The parallel case, given 
m Equation 6, and the series case given by 




k ' k 

solid gas 

•*—+( 1 -*)•*„ 


(9) 


are typically combined in one of two ways. One 
approach assumes the combined conduction thermal 
conductivity can be obtained by a superposition of the 
parallel and series arrangements 

k = A-k pM + {\-A)-k_ (10) 

where A is the fraction of heat transfer in parallel mode 
and ( 1 -A) is the fraction of heat transfer in series 
mode. 15 11 The other approach, used by Bhattaclirya, et. 
al., 10 assumes the combined conduction thermal 
conductivity is the square root of the sum of the squares 
of the parallel and series arrangements 

k = ^\k para J+{\-A)-{k ser J (11) 


The superposition of the parallel and series 
arrangements, Equation 10, with the solid thermal 
conductivity defined in Equation 7 will be used to 
model the combined conduction thermal conductivity in 
this study. Written in terms of the thermal conductivity 
of the gas and bulk material, the combined conduction 
thermal conductivity is given by 

k = A-[e k gas + (l-£)-F-k M ] 


+(i-4 


F ' k bulk ' k gas 




( 12 ) 


Depending on the values of A and F the forms of the 
combined conduction thermal conductivity given by 
Equations 6, 8, 9, and 1 0 are easily realizable. Setting 
/1=1 results in the parallel arrangement, while using 
A = 0 produces the series arrangement. The combined 
conduction thermal conductivity given by Equation 8 is 
obtained by letting A=l and F=(l-e) m l -Z). 


The thermal conductivity of the nitrogen gas 

k =^£L 

gas z 

was based on the gas conduction model used by 


(13) 


Daryabeigi, where k’ gas is the temperature dependent 


thermal conductivity of the gas at atmospheric 
pressure. 1 The term in the denominator, Z, is defined 
by 

Z = <S> + 'V-2?—^- Kn (14) 

a y + 1 Pr 


where a is the thermal accommodation coefficient, yis 
the specific heat ratio, Pr is the Prandtl number, and Kn 
is the Rnudsen number. The parameters <P and 0 
depend on the Rnudsen number. 0 = 1, 0= 0 for 
Rnudsen number less than 0.01 (continuum regime), 0 
= 1, 0= 1 for Rnudsen number between 0.01 and 10 
(transition regime), and 0=0, 0= 1 for Rnudsen 
number greater titan 10 (free-molecular regime). The 
Rnudsen number, Kn, is calculated h orn 
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(15) 


where A is the gas molecular mean free path and Jis die 
characteristic length. The mean free path is given by 


A = 


k b t 

n-Pd 1 

y g 


(16) 


where K B is the Boltzmann constant, T is die absolute 
temperature of the medium, P is the pressure, and d g is 
die gas collision diameter. For gas conduction in 
porous media, the characteristic lengdi is usually 
defined as die average linear pore size. 4 The pore size 
for metallic foam was assumed to have die same form 
as die pore size in fibrous insulation, giving a 
characteristic lengdi of 


S 


7Z D f 
~4 l-£ 


(17) 


where Df is the diameter of die strut. 1 The specific heat 
ratio, % for nitrogen gas is 7/5. The diermal 
accommodation coefficient, a, and the Prandtl number, 
Pr, for nitrogen gas were taken from Daryabeigi. 1 


G + - 


3 - [f 


dG_ 

dx 


= 4 a r; 


( 21 ) 


2 -e. 


2 J 


where T Ho , and T Co id are the temperatures at die 
boundary surfaces. The emittance of the septum plate, 

£ i, is 0.85 and die emittance of the water-cooled plate, 
£ 2 , is 0.92 as determined by Daryabeigi. 11 The 
assumptions used in diis formulation consist of: 
isotropic scattering, homogeneous and gray medium, 
diffuse emitting and reflecting surfaces. The 
assumption of isotropic scattering is not strictly valid, 
but provides simplification of the governing equations 
to yield an approximate solution. The extinction 
coefficient is determined from die specific extinction 
coefficient, e, by 

P = e-p (22) 


where p is die density of the foam. The specific 
extinction coefficient and die albedo of scattering are 
intrinsic properties of die material diat are independent 
of density and must be determined experimentally. 


Radiation 


Numerical Finite Volume Formulation 


Most of the work done on carbon foams and fibrous 
insulation has used die optically diick approximation 
for die radiation heat transfer. The optical diickness of 
the samples studied here was not known a priori; 
therefore, the use of the optically thick approximation 
could not be justified. The two- flux approximation, 
which is applicable over various ranges of optical 
thickness, was used to determine the amount of heat 
radiated through the void areas of the foam. ' The heat 
transferred by radiation, q ra d, is given by 


1 dG 
3 • [i dx 


(18) 


where fi is the extinction coefficient, the fraction of 
radiation lost to scattering and absorption per unit 
distance within the participating medium. The incident 
radiation per unit area, G, is determined from the 
incident radiation equation, 


G- 


1 d 2 G 
3 - fi 2 (l-w) dx 2 


= 4- a T 


(19) 


where ttris the albedo of scattering and <7 is the Stefan - 
Boltzmami constant. The incident radiation equation is 
subjected to the boundary condition at the septum plate 
of 


G — 


3 p 2 


.^ =4 .<rr;, 

^ dx 


2 -£ 


( 20 ) 


1 J 


and the boundary condition at the water-cooled plate of 


The analytical solution to the general governing 
equation of heat flow within the metallic foam is 
difficult to produce. Therefore, a numerical finite 
volume scheme was used on the transient heat flow 
problem. The thickness of the sample was separated 
into discrete volume elements and the conservation of 
energy principle was applied to each volume element. 
The change of the heat stored in each volume element 
plus the heat flowing out of each volume element was 
equal to the heat flowing into each volume element. 
The governing finite volume formulation, based on 
Equation 3, is given by 


p-c ■ Ax- 




At 


, T" -T" 1 G" -G" 

+ k. ■— — + - 1 — 


Ax 3 ■ /?; Ax 
_ k „ 7'" -7'/' | 1 G n ,-G" 


Ax 3 ■ /?", Ax 


(23) 


where the superscript denotes the tune step and the 
subscript denotes the spatial step. Boundary conditions 
applied were constant temperatures corresponding to 
the measured steady-state temperatures at the top and 
the bottom of the sample. A linearly varying 
temperature distribution through the thickness of the 
samples was selected as the initial condition. At each 
time step, the incident radiation equation. Equation 19, 
subjected to the boundary conditions of Equations 20 
and 21, was solved numerically using a finite difference 
technique to determine values of G based on the current 
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values of the temperatures through the thickness of the 
sample. The steady-state temperature distribution was 
found by iterating Equation 23 until the temperature at 
each volume element did not change with successive 
iterations. The effective thermal conductivity was 
calculated from Equation 1 using die temperature 
difference across die sample and die steady-state total 
heat flux obtained from die numerical solution. 


Parameter Estimation 


Discussion of Results 

The experimental data for die 14, 11, and 8 layer 
samples are presented and discussed. Then, die 
assumptions used in die development of die numerical 
model are explained. Using a subset of experimental 
data obtained from die 8 layer sample, parameters 
describing specific intrinsic properties of die foam are 
obtained. Results generated by the numerical model 
using die aforementioned intrinsic parameters are 
presented and compared to die experimental results. 


The unknown parameters, A and F, needed in die 
combined conduction diermal conductivity, Equation 
12, and die unknown parameters, e and to, needed in die 
two- flux approximation for die incident radiation 
equation were determined using die Levenberg- 
Marquardt mediod for nonlinear parameter estimation. 16 
The Levenberg-Marquardt mediod is an iterative 
procedure based on die minimization of the ordinary 
least squares norm given by 

S(P) = [K-Y(P)] r [K-Y(P)] (24) 


where P is a vector of N unknown parameters, K is die 
vector of M measured values of die effective thermal 
conductivity, and Y (P) is die vector of M predicted 
values of die effective thermal conductivity. The T 
superscript denotes the transpose. The iterative 
equations used in die Levenberg-Marquardt method of 
parameter estimation are 

A(P* +1 - Pj = b (25) 

where, 

A = (J t ) r J,+//A (26) 

is a NxN matrix and 

b = (J t ) r [K-Y(P t )] (27) 

is a N element vector. The subscript, k, denotes the 
iteration number. The sensitivity or Jacobian matrix 
coefficients are obtained from 


J ,.=- 


3T 

dP 


(28) 


where i = 1 to M and j = 1 to N. The damping 
parameter, //, along with the diagonal matrix, £2, damp 
oscillations and instabilities that arise out of the ill- 
conditioned nature of the problem. 


The solution to the system of linear algebraic equations, 
Equation 25, was used to produce a new set of 
parameters, P *+;. A negligible change in a successive 
set of calculated parameters served as the convergence 
criteria. 


Experimental Results 


Experimentally measured values of the effective 
diermal conductivity as a function of temperature 
difference across the sample for die 14, 1 1, and 8 layer 
samples are presented in Figures 4, 5, and 6, 
respectively. Data shown on die graphs are at nitrogen 
gas pressures of 10 4 , 1 .0, and 750 mm Hg. All diree 
samples behaved hi a similar fashion. At low pressure, 
where gas conduction was negligible, die primary 
modes of heat transfer were solid conduction and 
radiation. The magnitude of die effective thermal 
conductivity was at a minimum at low temperature 
differences. An increase in the temperature difference 
across the sample increased the contribution of 
radiation and the effective thermal conductivity 
increased as seen on the 10 1 mm Hg constant pressure 
curve. An increase in pressure increased die 
contribution of the gas conduction thus increasing die 
value of the effective thermal conductivity as evident in 
the 1.0 and 750 mm Hg constant pressure curves. 

Detailed error analysis was performed to determine die 
uncertainty due to the bias and random errors and 
spatial variations of the measured temperatures and heat 
fluxes. 1 7 The uncertainty hi the measured values of the 
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Figure 4. The effective thermal conductivity for the 
14 layer sample. 
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Figure 5. The effective thermal conductivity for the 
11 layer sample. 

effective thermal conductivity, depicted in Figures 4, 5, 
and 6 by the error bars, was within a 95% confidence 
level. At the lowest temperature differences, the error 
limits were much larger than the effective thermal 
conductivity. Neglecting those lowest temperature 
differences, the error limits for die 14 layer sample 
varied between 6 percent and 10 percent widi an 
average error limit of 8 percent. Similarly, die 1 1 layer 
sample had an average error limit of 4 percent and 
varied between 2 percent and 9 percent. The 8 layer 
sample error limits varied between 2 percent and 17 
percent with the average being 5 percent. 


not nearly as close as desired. The 14 layer and die 8 
layer data were significantly different from each other. 
The first 2 points of the 1 1 layer matched die 8 layer 
sample but die remaining points were more consistent 
with die 14 layer sample. Prolonged exposure of the 
layers near die septum plate to extreme temperatures 
diroughout die experimental investigation is expected to 
be the primary source of the deviations. The 8 layer 
sample was tested first dn oughout die entire 
temperature range dien used as a core for die 1 1 layer 
sample. The 1 1 layer sample was in turn used as a core 
for die 14 layer sample. 
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Figure 7. Comparison of the effective thermal 
conductivity of the three samples. 
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Figure 6. The effective thermal conductivity for the 
8 layer sample. 


A comparison of the effective thermal conductivity of 
the different samples is given in Figure 7. At the gas 
pressure of 10 4 and 1.0 mm Hg, the effective thermal 
conductivity of the three samples was nearly identical 
and well within the experimental uncertainty range, as 
expected. The effective thermal conductivity had the 
same value independent of the thickness of the sample. 
However, at the gas pressure of 750 mm Hg, the 
effective thermal conductivity of the three samples was 


Inspection of the final sample upon completion of the 
experimental investigation revealed the 6 layers closest 
to die septum plate were mechanically interlocked or 
stuck together. The interlocked layers would be 
expected to have a higher thermal conductivity than 
unlocked layers thus contributing to die higher 
measured values of the effective thermal conductivity 
for the 1 1 and 14 layers samples. 

Although the thickness of the layered samples was 
determined to within 0.03 nun. the distance between die 
septum plate and the water-cooled plate could only be 
determined to within 0.5 nun due to the surface 
variation of the plates. This uncertainty propagated 
dn ough the detailed error analysis and had a significant 
effect on the uncertainty of a single measured value. 
However, die dominant factor in the determination of 
die error limits depicted in the graphs was die 
uncertainty due to the spatial variation of the four 
measurement locations. The uncertainty due to die 
spatial variation was typically 2 to 3 orders of 
magnitude greater than the uncertainty of a single 
measured value. 
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Development of the Numerical Model 

For a selected set of boundary conditions (the 
temperatures at the top and the bottom of the sample), 
the steady-state temperature distribution and heat flux 
through the thickness of the samples were obtained 
from tlie numerical finite volume formulation of the 
analytical model of die heat transfer through die foam. 
Equation 23. Using die temperature difference across 
the sample defined by the boundary conditions and 
steady-state heat flux from die numerical solution, die 
effective thermal conductivity was calculated using 
Equation 1. 

The samples were constructed widi an integer number 
of foam layers. Assuming a uniform distribution of 
material tiiroughout die sample would simplify the 
analytical model; however, closer observation of the 
experimental data at a pressure of 10 4 and 750 mm Hg 
revealed some interesting trends diat could only be 
explained by incorporating die discrete nature of the 
layers into die analytical model. 

At 10 1 mm Hg, radiation and solid conduction were die 
main modes of heat transfer. Increasing die pressure to 
750 mm Hg essentially superimposed gas conduction to 
die data at 10 4 mm Hg. The difference between die 
data at the two pressures should be almost equal to die 
gas thermal conductivity over the temperature range; 
however, the difference was typically 5 times diat of the 
gas thermal conductivity. Furthermore, the effective 
thermal conductivities at 10 4 mm Hg were too low to 
incorporate the contribution of the solid conduction, hi 
other words, a solid conduction efficiency factor, used 
in Equation 12, could be found to account for the solid 
conduction at 10 1 mm Hg but the same efficiency 
factor produced a gross underestimation at 750 mm Hg. 
If perfect thermal contact between the different layers 
in each sample and between the outermost layers and 
bounding septum and water-cooled plates is assumed, 
the solid conduction contribution should produce 
effective thermal conductivity values much higher than 
what was measured at low pressure. 

The only physical explanation for the observed 
phenomena is that there was a discontinuity in solid 
conduction between the sample layers and between the 
outermost layers and the bounding plates (septum and 
water-cooled). The only conduction mechanism 
between the sample layers and between the outermost 
layers and the bounding plates was gas conduction. 
Thus, at 1 0 1 mm Hg, since gas conduction was 
negligible, there was little appreciable transfer of heat 
between the sample layers by conduction, and the net 
heat flux through the sample was due mostly to 
radiation. As gas pressure increased, heat was 


conducted between the sample layers and between the 
outermost layers and the bounding plates by gas 
conduction. 

This discontinuity in solid conduction was modeled by 
placing a small gap between the sample layers and 
between the outermost layers and the bounding plates. 

A gap spacing of 0.01 mm was used. The gap contained 
only gas, thus at low pressure the primary mode of heat 
transfer across the gap was radiation. The gap was 
incorporated into the numerical finite volume 
formulation by the use of a non-uniform node spacing 
across the thickness of the samples. The location of the 
nodes was selected so a gap corresponded to an integral 
number of nodes. The conduction thermal conductivity 
through the volume elements, defined by the gap nodes, 
was due to gas conduction only, Equation 13. Within a 
single layer, it was assumed the combined conduction 
thermal conductivity of Equation 12 could be used to 
adequately describe the behavior of the metallic foam. 

The solution to the incident radiation equation was 
applied across the thickness of the sample based on the 
same non-uniform node spacing used above. However, 
the density distribution was assumed uniform 
throughout the sample for the radiation calculations. 

The extinction coefficient, albedo of scattering, and the 
solid conduction efficiency factor are parameters 
intrinsic to the material under investigation; however, 
these parameters were not known. Therefore, a subset 
of measured values was used to predict these 
parameters based on the Levenberg-Marquardt method 
for nonlinear parameter estimation. 

The subset of measurements chosen for use in the 
parameter estimation method was selected from the 8 
layer sample at 10 1 and 750 mm Hg. At high vacuum, 
10 4 mm Hg, the gas was within the free-molecular 
regime and conduction through the gas within the void 
of the foam had a minimal effect, hi this regime the 
contribution of the combined conduction thermal 
conductivity was expected to have an insignificant 
influence on the effective thermal conductivity. For 
this reason, measurements at high vacuum were used to 
estimate the radiation dependent parameters, e and ox 
At high pressure, 750 mm Hg, the gas was within the 
continuum regime defined by the temperature 
dependent thermal conductivity of the gas at 
atmospheric pressure, hi this regime the contribution of 
the combined conduction thermal conductivity was 
expected to have a significant contribution to the 
effective thermal conductivity. The measurements at 
high pressure were used to estimate the weighting 
factors for modeling the combined conduction thermal 
conductivity. 
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Numerical Results 

The numerical finite volume formulation produced die 
same values of the effective thermal conductivities for 
die 14, 11, and 8 layer samples. Since die numerical 
results are independent of the number of layers (dius 
die diickness of the sample) only die 8 layer sample 
results are presented. 

Using the combined conduction tiiermal conductivity 
given by Equation 12 in die numerical model produce 
die curves shown in Figures 8 and 9. The total heat 
flux dirough die sample is shown in Figure 8 for die 
constant temperature difference across die sample of 
542 K and nitrogen gas pressure of 750 mm Hg. The 
total heat flux was constant as expected for a steady- 
state condition. The curves were generated for die 8 
layer sample using 5 nodes across each layer. Also 
shown on die graph are die components of the total heat 
flux, die heat flux due to combined conduction and die 
heat flux due to radiation. The discontinuous jumps in 
heat flux, observable at every fifth node, were due to 
the presence of gaps between layers, hi die figure, die 
septum plate was located at node position 0.0 mm and 
the water-cooled plate was located at 13.69 mm. The 
dominant mode of heat transfer across die diickness of 
the sample was combined conduction. Radiation 
became more significant in die proximity of the septum 
plate. 



0 5 10 15 

Node Position (mm) 


Figure 8. Spatial variation of the total, conductive, 
and radiative fluxes for the 8 layer sample at 750 
mm Hg. 

The measured and numerically predicted effective 
thermal conductivities of the 8 layer sample are shown 
in Figure 9 for the nitrogen gas pressures of 10 4 and 
750 mm Hg. The numerical results ate presented by the 
solid curves. The radiation dependent parameters were 
found to be e = 10.23 — 1.77x10 3 T and a>= 0.8 where 
T is the absolute temperature. For A = 1, analogous to 
the parallel arrangement, the solid conduction 


efficiency factor was determined to be F = 0.061. The 
identical curve was also produced for A = 0.802 and F 
= 0.080, die parallel-series arrangement. As seen from 
die figure, die numerical results had excellent 
correlation to experimental results at low pressure but 
failed to capture die behavior at high pressure. 



0 200 400 600 800 

AT(K) 

Figure 9. Comparison of measured and predicted 
effective thermal conductivities using Equation 12 
for the combined conduction thermal conductivity. 

The difference in effective thermal conductivity 
between die high pressure and low pressure curves 
produced by the numerical model was approximately 
constant while the experimental data had an increase in 
the difference as the temperature difference across the 
sample increased. The combined conduction thermal 
conductivity model given by Equation 12 assumed a 
superposition of the solid conduction and the gas 
conduction modes. There must be a coupling or 
interaction between the solid and gas conduction diat is 
not accounted for in Equation 12. To correct for die 
discrepancy a coupling term was introduced into die 
numerical model. The coupling thermal conductivity 
used, k coupling , is given by 

k =a\k gas -F- k bult ) 2 (29) 

where a is a coupling weighting factor. Since die 
thermal conductivity of gas is extremely sensitive to 
both pressure and temperature, the coupling diermal 
conductivity is highly dependent on pressure and 
temperature. The resulting combined conduction 
tiiermal conductivity, assuming a parallel arrangement 
(A = 1) in Equation 12, is 

k = £-k sas + (l-£)-F -k bult + a-{k gas F k bulk Y (30) 

When Equation 30 was used for the combined 
conduction thermal conductivity in the numerical 
model, the curves shown in Figure 10 were produced. 
For A = 1. die parameter estimation procedure 
generated die parameters of F = 6.85x10 and a = 
389.0. The radiation dependent parameters were 
determined to be e = 9.85 — 2.63x10° T and a>= 0.993. 
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The radiation dependent parameters were slightly 
affected by the addition of the coupling term since a 
change in the solid conduction efficiency factor was 
required. As seen from the figure, the addition of the 
coupling term placed the numerical results within 
experimental uncertainties for both the 10 and 750 
mm Hg regions. 



Figure 10. Comparison of measured and predicted 
effective thermal conductivities using Equation 30 
for the combined conduction thermal conductivity. 

The gas conduction model given by Equations 13 
through 17 was developed for porous media with 
uniform density fibers or struts. As seen in Figure 2, 
however, the strut of the metallic nickel foam is hollow. 
Therefore, the thickness of the strut wall was used in 
the place of the diameter of the strut hi Equation 17. 

The thickness of the strut wall was approximately 0.014 
mm. The calculated values of the effective thermal 
conductivity for the 8 layer sample with a temperature 
difference of 500 K across the thickness of the foam 
over the pressure range of 10 1 to 750 mm Hg are 
shown in Figure 1 1. Although the predicted values do 
not correspond within experimental error limits within 
the region defined by the pressure of 0.001 to 



0.0001 0.001 0.01 0.1 1 10 100 1000 

Pressure (mm Hg) 


Figure 11. Variation of effective thermal 
conductivity with pressure. 


0. 1 mm Hg, the predicted values do show the general 
trend throughout the rarefied, transition, and continuum 
regions. Neglecting the region between 0.001 and 0. 1 
mm Hg, an average 9 percent difference was obtained 
between the calculated and measured values. The 
largest percent difference occurred at low pressure, with 
almost 22 percent difference, and the smallest percent 
difference occurred at the intermediate pressures, with 
less than 1 percent difference. 

Conclusion 

The effective thermal conductivity of metallic nickel 
foam was measured for a temperature difference range 
from 100 to 1000 K and a pressure range of 10 4 to 750 
mm Hg. A numerical model was developed to predict 
the behavior of the effective thermal conductivity at 
various temperatures and pressures. Using a small 
subset of experimental data, parameters related to 
intrinsic properties of the foam were determined. The 
calculated values of the effective thermal conductivity 
of the metallic foam produced by the numerical finite 
volume formulation using the predicted intrinsic 
parameters were compared to die measured values not 
used in the parameter estimation metiiod to validate die 
numerical models. Correlation of die numerical results 
to die experimental values required the introduction of a 
conduction coupling term to the gas/solid conduction 
model. Calculated values corresponded to within an 
average of 9 percent to die experimental values. The 
model was consistent with experimental results 
throughout die environmental conditions under 
examination. 


Acknowledgements 

Research done by the author is a portion of what will be 
included in a Master of Science Thesis widi The 
George Washington University. The authors would like 
to thank Mr. Jeffery R. Knutson from NASA Langley 
Research Center for providing temperature and pressure 
control and his assistance in conducting the tests, and 
Mr. James M. Baughman horn AS&M at NASA 
Langley Research Center for producing the 
photomicrograph images. 


References 

1. Daryabeigi, K., “Design of High Temperature Multi- 
layer Insulation for Reusable Launch Vehicles”. Ph.D. 
dissertation, LTniversity of Virginia, May 2000. 


10 

American histitute of Aeronautics and Astronautics 





2. Fuganti, A., Lorenzi, L., Hanssen, A.G., and 
Langseth M., “Aluminum Foams for Automotive 
Applications,” Metal Foams and Porous Metal 
Structures, ed. by Banhart, J., Ashby, M.F., and Fleck, 
N.A., Verlag MIT Publishing, 1999. 

3. Ettel, V.A., “New Inco Powders and Foams for 
Nickel Batteries,” NiCad 98, Prague, Sept. 21-22, 1998. 

4. Kaviany, M., Principles of Heat Transfer in Porous 
Media, 2 nd ed., Springer- Verlag New York, hie, 1995. 

5. Lee, S.-C., and Cunnington, G.R., “ Conduction and 
Radiation Heat Transfer in High Porosity Fiber 
Thermal Insulation,” Journal of Thermophysics and 
Heat Transfer, Vol. 14, No. 2, April-June 2000, pp. 121- 
136. 

6. Glicksman, L. R., “Ch. 5, Heat Transfer in Foams,” 
Low Density Cellular Plastics, ed. by Hilyard, N.C., 
and Cunningham, A., Chapman and Hall, 1994. 

7. Gibson, L.J., and Ashby, M.F., Cellular Solids, 
Structures and Properties - Second Edition, Cambridge 
University Press Ltd., 1997. 

8. Baillis, D., Raynaud, M., and Sacadura, J.F., 
“Determination of Spectral Radiative Properties of 
Open Cell Foam: Model Validation”, Journal of 
Thermophysics and Heat Transfer, 14, No. 2, April- 
June 2000. 

9. Calmidi, V. V. and Mahajan, R.L., “The Effective 
Thermal Conductivity of High Porosity Fibrous Metal 
Foams”, ASME Journal of Heat Transfer, 122, May 
1999, pp. 466-471. 

10. Bhattacharya, A., Calmidi, V.V., Mahajan, R.L., 
“An Analytical-Experimental Study for the 
determination of the effective thermal conductivity of 
high porosity fibrous foams” , Application of Porous 
Media Methods for Engineered Materials, AMD- Vol. 
233, ASME 1999, pp. 13-20. 

11. Daryabeigi, K., “Analysis and Testing of High 
Temperature Fibrous Insulation for Reusable Launch 
Vehicles”, 37th AIAA Aerospace Sciences Meeting and 
Exhibit, AIAA 99-1044, January 11-14, 1999. 

12. Tannehill, J.C., Anderson, D.A., and Pletcher, R.H., 
Computational Fluid Mechanics and Heat Transfer, 
Second Edition, Tayler & Francis, 1997, pp. 72. 

13. Bauer, T.H., “A General Analytical Approach 
toward the Thermal Conductivity of Porous Media,” 

Int. J. Heat Mass Transfer, Vol. 36, No. 17, 1993, pp. 
4181-4191. 

14. Kreith, F. and Black, W.Z., Basic Heat Transfer, 
Harper & Row, Publishers, 1980, pp. 509. 


15. Pawel, R.E., McElroy, D.L., Weaver, F.J., and 
Graves, R.S., “High Temperature Thermal Conductivity 
of a Fibrous Alumina Ceramic,” Thermal Conductivity, 
19, ed. by David Yarbrough, Plenum Press, New York, 
1985, pp. 301-313. 

16. M. Necati Ozisik, Helcio R. B. Orlande, Inverse 
Heat Transfer, Taylor & Francis, 2000, pp. 37. 

17. Coleman, H.W., Steele, W.G., Jr., Experimentation 
and Uncertainty Analysis for Engineers, Joint Wiley & 
Sons, hie, 1989. 


11 

American histitute of Aeronautics and Astronautics 



